High energy protons from PKS 1333-33 
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In this letter we give an account of the possible acceleration of protons in the outer radio lobes of 
the active galaxy PKS 1333-33. We also make estimates of the arrival energy spectrum. 
PACS number(s): 96.40 - 98.70.S - 95.85.R - 13.85. T 



Radio galaxies considered as sources of cosmic rays 
(CR) beyond the Greisen-Zatsepin-Kuz'min cutoff |] 
must be quite close to our own galaxy (z < 0.03). There 
are just a few objects within this range and, consequently, 
the expected ultra high energy CR distribution on the 
sky should be highly anisotropic. Assuming that the in- 
tergalactic magnetic field is near 1 nG, it should be ex- 
pected an excess of CR detections at energies larger than 
5 x 10 19 eV within regions of angular radius 9 < 15° 
and centered at the positions of the nearest active galax- 
ies According to this picture, the southern CR-sky 
should be dominated by three outstanding sources: Cen 
A (the nearest radio galaxy) which would provide the 
most energetic particles detectable on Earth B , Pictor A 
(a strong source with a flat radio spectrum) which would 
contribute with the larger CR flux [§), and PKS 1333-33 
which might be a source of events similar to those re- 
cently detected in the Northern Hemisphere. In addition 
to these sources, there are other two southern candidates, 
Fornax A (z = 0.057) and PKS 2152-69 (z = 0.027), 
which could provide contributions to the CR flux above 
the cutoff. An approximate theoretical picture of the 
up to now unexplored ultra high energy CR southern 
sky is consequently at our disposal. The present letter 
is devoted to CR production in the outer region of the 
nearby southern active galaxy PKS 1333-33 and the spec- 
tral modifications arisen from their propagation through 
intergalactic space. 

The radio galaxy PKS 1333-33 is made up of a core, 
two symmetric jets, and two extended radio lobes. Fig. 
1 shows the entire source at 20 cm according to the VLA 
map obtained by Killeen et al. ||. The core has been 
identified with the El galaxy IC 4296, which has a red- 
shift z=0.013 The distance to the source is 35.2 h^ 1 
Mpc, (h = iJo/100 km -1 s Mpc, hereafter we shall adopt 
h = 0.65 §). The large-scale structure of PKS 1333-33 
has been studied in detail at 1.3, 2, 6, and 20 cm with an 
angular resolution of 3.2" ||. The jets are slightly bent, 
presumably as a consequence of the motion of the core 
with respect to the intergalactic medium. The total flux 
density of the source at 20 cm is ~ 14.5 Jy. The inte- 
grated radio luminosity is ~ 2.5 x 10 41 erg s _1 assuming 
a spectral index of a = 0.6 and frequency cutoffs at 10 7 
and 10 11 Hz. Actually, the value a = 0.6 is correct only 
for the extended radio lobes. The spectral index steep- 



ens in the jets, reaching a = 1, while the core has flatter 
values: a ~ 0.3. 

The hot spots of extended radio sources are regions of 
strong synchrotron emission These regions are pro- 
duced when the bulk kinetic energy of the jets ejected by 
a central active source (supermassive black hole plus ac- 
cretion disk) is reconverted into relativistic particles and 
turbulent fields at a "working surface" in the head of the 
jets ||. The speed i>h with which the head of a jet ad- 
vances into the intergalactic medium of particle density 
n e can be obtained by balancing the momentum flux in 
the jet against the momentum flux of the surrounding 
medium. Measured in the frame comoving with the ad- 
vancing head, Vh w Uj [1 + ( J ie/ rl j) 1 ^ 2 ]~ 1 , where n,j and i>j 
are the particle density and the velocity of the jet flow, 
respectively. Clearly, v s > «h for n e > rij , in such a way 
that the jet will decelerate. The result is the formation 
of a strong collisionless shock, which is responsible for 
particle reacceleration and magnetic field amplification. 
The acceleration of particles up to ultrarelativistic ener- 
gies in the hot spots is the result of repeated scattering 
back and forth across the shock front. The particle de- 
flection in this mechanism is produced by Alfven waves 
in the turbulent magnetic field. Biermann and Strittmat- 
tcr [[l0| have studied this process assuming that the en- 
ergy density per unit of wave number of MHD turbulence 
is of Kolmogorov type. According to their calculations, 
the highest energy of protons injected in the intergalactic 
medium from the hot spot can be obtained by balancing 
the gains and losses in the diffusive shock acceleration 
process: 

£ max = 7.8 x 10 5 u 3 ^ R- 1 ' 2 (1) 
x £Q 5/4 (l + Aa)~ 3 / 4 EeV 

where /3 ot stands for the jet velocity in units of c, u is 
the ratio of turbulent to ambient magnetic energy den- 
sity in the hot spot (of radius i?hs measured in kpc), -B_5 
is the magnetic field in units of 10~ 5 G, a is the ratio of 
photon to magnetic energy density, and A gives a rela- 
tive strength of jp interactions against the synchrotron 
emission. 

An interesting feature of PKS 1333-33 is an intense re- 
gion of synchrotron emission localized at the outer edge 
of the eastern lobe. This region can be considered as a 
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FIG. 1. VLA image of the entire PKS 1333-33 source ob- 
tained at a wavelength of 20 cm by Killeen et al. 



The radius of the acceleration region can be directly 
measured by means of a Gaussian fitting from the de- 
tailed VLA maps obtained by Killeen et al. @J, result- 
ing i?h s ~ 2.5 hr 1 kpc. The velocity of the jet is not 
well established. If the source is ~ 10 7 yr old, a velocity 
/3jet ~ 0.2 can be estimated from an analysis of the energy 
budget |ll| . The value of A in Eg. (§) has been computed 
by Biermann and Strittmatter [ |10| almost independently 
of the source parameters. They obtain A ps 200. A value 
a ~ 0.01 seems to be reasonable for a source with the 
luminosity of PKS 1333-33 [ |l0] |. Taking the above con- 
siderations into account (with a typical value for the total 
magnetic field in the spot of £>_5 ~ 10 [Q), we obtain 
from Eq. (0) the maximun injection energy for protons, 
£max « 6 x 10 20 eV. 

The diffusive shock acceleration process in the working 
surface of the eastern lobe leads to a power law particle 
spectrum [13], 



Q(E) = k E- T (E <E< £ max ) 



(5) 



"working surface" formed by the deceleration of the jet. 
This interpretation is supported by VLA polarimetric ob- 
servations, which show a change in the field orientation 
from parallel to perpendicular to the jet axis ||. This 
change is probably due to the rearrangement of the field 
lines in the post-shock region. The degree of linear polar- 
ization in the eastern lobe is in the range 20%-40%. The 
synchrotron parameters estimated for this region are: the 
minimum energy density (e m i n = 2.1 x 10~ 13 erg cm -3 ), 
the minimum magnetic field (B min = 1.5 x 10 -6 G), the 
minimum total pressure (P m - m = LI x 10~ 13 dyne cm -2 ), 
and the degree of linear polarization (m = 20-40%). 

The leptonic component of the CRs will produce syn- 
chrotron emission with a spectrum given by S u oc v~ a , 
where a = (F - l)/2. Since a = 0.6, we get F = 2.2. In 
what follows, we shall assume that electrons and protons 
in the source obey the same power law energy spectrum 
oc£- r . 

The degree of linear polarization expected for the syn- 
chrotron radiation when the magnetic field is homoge- 
neous is 



m (r) 



3L + 3 
3F + 7 



71% 



(2) 



However, the observed degree of polarization has a mean 
value of m(r) ~ 30%. This fact can be explained by the 
presence of a turbulent component B r in the field, in such 
a way that 



m(r) = TOo(r) 



(3) 



where Bq stands for the homogeneous field. From Eq. 
(§) and (§) we get 



B r ps 1.2 B 
and consequently u — -B 2 /_B 2 otal ps 0.6. 



(4) 



where k — (no/P — 1) E^ ~ r for E max 3> Eq, and no is 
the particle density in the source. Thus, we have the 
following proton injection spectrum 



Q(E) cx E 



-2.2 



E < E < 6 x 10 20 eV 



(6) 



Let us now estimate the losses that the high energy 
protons suffer during their journey through the inter- 
galactic medium, and the consequent modification of the 
injection spectrum. This energy degradation originates 
in the interaction of the protons with the very low energy 
photons of the cosmic microwave background (CMB). For 
nearby sources (distances < 100 Mpc) photomeson pro- 
duction is the dominant interaction at E > 10 19 eV (if 
one considers the contribution coming from the tail of the 
Planck distribution) jl4j . The propagation of the protons 
through the relic gas of photons may be modelled by a ki- 
netic transport equation. Although it is clear that the in- 
teractions of the ultra high energy protons with the CMB 
photons lead to a step-by-step energy loss (which should 
be included in the transport equation as a collision in- 
tegral) , hereafter we shall use the continuous energy loss 
approximation assuming straight line propagation which 
is expected to be in agreement with the physics of the 
problem at hand. 

Using the formalism presented in |fj|] it is straightfor- 
ward to compute the main characteristic of the evolved 
spectrum. It is well known that the conservation of the 
differential particle number implies that, 



dN _ d(b(E)N) 
~dt ~ dE 



D V 2 N + Q, 



(7) 



where, in the first term on the right, b(E) is the 
mean rate at which particles lose energy (ps 3.66 x 
10 -8 E exp{287/£;} EeV/yr @); the second term, the 
diffusion in the CMB, is found to be extremely small due 
to the low density of relic photons and is neglected in 
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FIG. 2. Modification factor of PKS 1333-33. 

the following. The third term corresponds to the particle 
injection rate into the intergalactic medium given by Eq. 
(^). The solution of Eq. (0) reads, 



,20 



N(E,t) = kE„ 



b(E g ) 
b(E) 



(8) 



with the injection energy E g fixed by the constraint, 



3.66 x HT^yr- 1 - Ei(287/£) 



Ei(287/£; g ) = (9) 



being Ei(x) the exponential-integral function. The evo- 
lution of the injection spectrum can be then conveniently 
represented by the modification factor r), defined as the 
ratio between the modified and the unmodified spectrum. 
In Fig. 2 we have plotted the relation obtained for rj in 
the case of PKS 1333-33. The energy loss by photomeson 
production creates the expected cutoff, and the result- 
ing ultra high energy CRs (protons and neutrons) pile 
up just below the threshold energy of photopion produc- 
tion, forming a bump. However, the energies at which 
the cutoff and the bump appear seem to indicate that 
CRs above 100 EeV might be expected from the location 
in the sky of PKS 1333-33 (I « 313.7°, b w 27.7°). 

The question whether this scenario is the correct one 
or not should be answered in a few years by the Pierre 
Auger Southern Observatory [15| (fluorescence detector 
plus ground array) as well as by the future eyes of the 
OWL fll that will deeply watch into the CR-sky. 
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